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ABSTRACT 


The East Asian winter monsoon cold surge is examined using Navy tropical Global 
Band Analysis data for January-March of 1975-1988. Multiple-set Canonicle Correlation 
(MCC) analysis, single-point correlation analalysis, and composite analysis are used to 
determine the the cold surge structure over a broad domain that includes the western 
Pacific and eastern Indian ocean. 

The results show that the structure of the cold surge is manifested in MCC modes 
1 and 2 of surface v calculated over both the East and South China Seas. Single-point 
correlation and composite analyses indicate that the surge can be viewed as a single mode 
over a rather broad region, with an average period of 8 days, that includes simultaneous 
eastward propagating longwaves across the East China Sea and northwestern Pacific, and 
equatorward propagating surge waves across the South China Sea. Composite analyses 
during strong events also reveal cross-equatorial surge influence present at low levels and 
increases in divergence over the maritime continent at 200 hPa. 

MCC modes 3 and 4 calculated over a small spatial domain that includes the East 
China Sea describe 3-day eastward propagating waves with wavelengths on the order of 
1500-2000 km. Associated with these shorter waves are modest surge waves that also 


propagate equatorward across the South China Sea. 
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I. INTRODUCTION 


The piirpose of this paper is to examine the statistical characteristics of some 
tropical and subtropical synoptic-scale phenomena occurring during the East Asian winter 
monsoon based on a 14-year dataset. The statistical analyses are carried out using the 
Multiple-set Canonical Correlation Analysis (MCCA) technique that was recently 
developed by Chen et al. (1994). A brief review of the history of synoptic tropical 
motions is provided in section A. In section B, two significant aspects of the East Asian 
winter monsoon, broad-scale cold surges that penetrate into the tropics, and cyclogenesis 
in the East China Sea, are reviewed. Chapter II contains a brief description of the data 
sets and an overview of the analysis techniques used including MCCA. Chapter III 
provides a discussion of results obtained when MCCA is applied to two separate domains: 
the East China Sea and the South China Sea. In Chapter IV, composite analysis and 
single-point correlation analysis are performed based on the structure of the leading 
MCCA modes for the two regions. Chapter V provides a summary and conclusions. 

A. SYNOPTIC DISTURBANCES IN THE TROPICS 

Although the existence of wave disturbances in the tropics has been recognized 
since the 1930’s, detailed study of tropical synoptic processes was not possible until upper 
air soundings became widely available after World War II. Prior studies were done using 
only surface data (e.g., Piersig 1936; Regula 1936; Dunn 1940). Upper air data from the 
Caribbean area was used by Riehl (1945) in his "waves in the easterlies" model which 
describes westward moving, wave-like oscillations in the basic easterly current of the 
lower troposphere. Even with upper air soimdings, the study of tropical dynamics was, 
and still is, difficult due to the limited amount of data. 

In the late 1960’s, spectral analysis was beginning to be used to pull information 
from rawinsonde data collected at widely dispersed observing stations. Yanai et al. (1968) 
applied spectral analysis to rawinsonde data over the equatorial Pacific and found 
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westward propagating waves with eastward vertical tilt having periods of 4 to 5 days. 
These waves with wavelengths on the order of 6000 to 10000 km were subsequently 
identified with mixed Rossby-gravity waves in the equatorial stratosphere (Maruyama 
1967, 1968). Wallace and Chang (1969) and Chang et al (1970) also foimd waves in the 
equatorial western Pacific with a dominant 4 to 5 day period. These waves in the lower 
troposphere have a zonal wavelength of 20(X) to 4000 km and a vertical tilt that 
systematically changes from eastward in the central Pacific to slightly westward in the 
extreme western Pacific. This east-west variation in vertical structure was also seen later 
when comparing years with warm and cold sea surface temperatures over the central and 
eastern Pacific. Warm years show a more equivalent barotropic structure while cold years 
show a more easterly tilt with height (Chang and Miller 1978). 

In 1979, results of the First GARP Global Experiment (FGGE) fueled further 
investigations into the behavior of tropical disturbances over areas where there were 
previously little or no station data available. Nitta et al. (1985) and Nitta and Takayabu 
(1985) conducted a global survey of the tropical transients using FGGE observations for 
July and August. Nitta and Takayabu (1985) and Tai and Ogura (1987) also used the 
FGGE analyses to show the existence of west-northwestward propagating disturbances 
in the eastern Pacific. 

Since 1979, global datasets have been analyzed and archived at various operational 
centers around the world. The availability of these global datasets greatly enhances the 
ability to conduct empirical tropical studies. Studies are no longer restricted to sparse 
station data or observations taken within limited-area networks during intensive field 
experiments. Of course, niunerical model analysis data, while having global coverage, 
may not always be reliable, especially in the tropical data-sparse areas. Nevertheless, they 
have proven to be very useful. Reed et al. (1988) studied the synoptic evolution of 
African wave disturbances appearing in the European Center for Medium Range Weather 
Forecasts (ECMWF) analyses for the summer of 1985. They found that independent 
satellite data confirmed the African wave features depicted by ECMWF analyses over 
data-sparse regions. 
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More recently, Liebmann and Hendon (1990) studied the 1980-1987 ECMWF 
analyses of the near equatorial wind field during the September-December season. Their 
focus was primarily on the behavior of synoptic-scale disturbances occurring within 10 
degrees north and south of the equator. They foimd westward propagating and equatorial 
trapped features in their 850 hPa meridional wind field with periods less than 6 days. 
The features looked somewhat like mixed Rossby-gravity waves. Lau and Lau (1990) 
used twice daily ECMWF global gridded analyses for the northern summers of 1980- 
1987. They examined lag-correlation and regression statistics of tropical fluctuations with 
synoptic time scales. Their results confirmed that a substantial amount of synoptic-scale 
variability in the tropics during northern summer is associated with propagating wave-like 
disturbances that remain coherent for several days. 

In addition to the synoptic-period disturbances discussed above, intraseasonal 
oscillations have also been studied and found to have significant influence on tropical 
weather. The most weU-known of these is the Madden-Julian Oscillation (MJO, Madden 
and Julian 1971). This is an eastward propagating, planetary-scale wave with zonal 
wavenumber one and dominant oscillations in the east-west (u) wind component. 
Observational studies of the MJO have shown that it plays a key role in governing 
tropical convection, particularly in the Indian monsoon region (e.g., Lau and Chan 1983). 
Salby and Hendon (1994) investigated the spectral character of tropical convection in an 
11 year record of outgoing longwave radiation fi-om the Advanced Very High Resolution 
Radiometer. They found that the space-time spectrum of convection possesses a broad 
peak at wavenumbers 1-3 and eastward periods of 35-95 days. This broad signal is 
associated with a large-scale anomaly that moves across and modulates the "climatological 
convection" over the equatorial Indian Ocean and western Pacific. Cardenas Amores 
(1994) used Navy Operational Global Analysis and Prediction System (NOGAPS) data 
for the summers 1989-1991 to study intraseasonal oscillations other than the MJO. Using 
MCCA, he foimd a near 20 day oscillation with a zonal half-wavelength around 3000- 
4000 km. Chen and Chen (1993) found a similar 10-20 day intraseasonal monsoon mode 
when they studied the 1979 FGGE data and daily rainfall estimates over India. 
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Most of these previous observational studies of the tropical synoptic scale motions 
focused on the northern summer period. In this study the focus will be on the northern 
winter. 


B. THE EAST ASIAN WINTER MONSOON 

1. Cold Surges 

Perhaps the most well known event that occurs dming the East Asian monsoon 
is the cold surge. According to Boyle and Chen (1987), a cold surge starts when 5(X) hPa 
northwesterly flow over Lake Baikal sets up the larger scale flow upon which short waves 
propagate toward the mean longwave trough that typically resides along the East Asian 
coast. The shortwaves move into phase with and amplify the longwave which produces 
a circulation that pushes cold air southward from eastern Siberia. At the siuface, an 
anticyclone develops over China as a result of the low-level cold air advection and upper- 
level vorticity advection. The result is an outbreak of very cold air that flows eastward 
off the China coast and southward into the South China Sea that can have a significant 
impact on regional weather (Ramage 1975). The cold surge generally follows isentropic 
surfaces so that by the time it reaches the South China Sea, the surge is relatively shallow 
and its signal is usually not detected above 850 hPa. Although there exist many 
definitions of surge criteria, general consensus is that the event is marked by rapidly 
dropping surface temperatures over southern China accompanied by an increase in the 
northeasterly wind both off the East Asian coast and in the South China Sea (Boyle and 
Chen 1987). 

The cold surge events are a good example of the interaction that takes place 
between the midlatimdes and the tropics. Chang and Lau (1982) analyzed 2(X) hPa winds 
of four winters and found that during surge periods, the midlatitude and tropical 
circulation components vary in a coherent way. Several investigators have noticed that the 
cold surge is often followed by increased tropical convection over the maritime continent 
(e.g., Lau and Chang 1987). Smdies of data taken during the Winter Monsoon Experiment 
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(WMONEX) for the winters 1978-1979 have shown that the degree of convective activity 
in the southern South China Sea is modulated significantly by diurnal forcing of land/sea 
breeze circulations (Houze et al. 1981; Johnson and Priegnitz 1981). Chu and Park (1984) 
studied a cold singe event from the WMONEX and found a rapid southward mass 
transport in the lower troposphere from the midlatitudes to the tropics. Chang et al. 
(1979) studied the structures of cyclonic vortices near Borneo and found that the cold 
surge actually deepens the vortex probably due to increased organization of cumulus 
convection. Chang and Chen (1992) conducted a statistical study of 103 surge cases 
during 14 winter seasons. They found a relationship between the occurrence of cold 
surges over the South China Sea and the variation of the 200 hPa and smface tropical 
divergence from the western Pacific to the eastern Indian Ocean. 

In this study we use the MCCA method to isolate the horizontal structure of the 
cold surge through surface meridional wind analysis. The time variation of the three- 
dimensional structure is then deduced from single-point correlation and composite 
analyses. 

2. Cyclones and Cyclogenesis in the East China Sea 

The nearly stationary major longwave in the northern hemisphere is located along 
the eastern coast of Asia. Manabe and Terpstra (1974) showed that the longwave pattern 
observed is probably due to topographic effects. Although the presence of the longwave 
has little effect on the growth rate and phase speeds of the fastest-growing linear 
disturbances, its location plays a key role in producing regions favorable to cyclone 
development downstream (Frederikson 1979). Observation of the oceanic regions east of 
the Asian continent, including the Sea of Japan and the East China Sea, certainly bears 
this out. Most of the Pacific ocean storms begin in these regions. Weak shallow 
disturbances that move off of southern mainland China often deepen and may deepen 
rapidly if an upper level trough reaches the coast coincidentally. 

Nitta and Yamamoto (1973) analyzed the cyclogenesis frequency over East Asia 
(20°N to 45°N, 105°E to 150°E) for the period 1966 to 1970 and identified two types of 
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cyclonic disturbances. The first type is the typical baroclinic wave in which growth is due 
to barochnic instability triggered by an upper-level trough in the westerlies with 
wavelengths on the order of 3000-4000 km. The second type is an intermediate-scale 
wave that has a horizontal wavelength of 1000 to 2000 km and is mostly confined to the 
lower half of the troposphere. At the initial stage of formation, the disturbance does not 
couple with any trough in the upper troposphere. Nitta and Yamamoto (1974) analyzed 
the structure of these intermediate-scale disturbances and found that the most frequently 
observed type in the winter has the warm air to the east side of the cyclone and cold au¬ 
to the west in the lower troposphere, with cold air occupying the midtroposphere above 
the cyclone center. Boyle and Oien (1987) summarized the characteristics of the 
mtermediate-scale disturbances originally given in a report for the Air Mass 
Transformation Experiment (AMTEX) (GARP 1981) as foUows: 

1. The characteristic scale of the disturbance is 1000 to 2000 km in the zonal 
direction. 

2. The characteristic time scale is short when compared with that of the typical 
baroclinic wave. 

3. The disturbances become active in a moist lower troposphere under the 
condition of a less stable stratification. 

4. These disturbances do not appear to be directly associated with an upper 
tropospheric trough, and their kinetic energy is confined to the lower portion of 
the troposphere. 

5. Accordmg to the radar observations, cumulus activity seems to be closely 
associated with the generation and development of these dismrbances. 


Boyle and Chen (1987) also note that the cyclones apparently require heat and 
moisture fluxes that are found in the East China Sea in order to develop and be 
mamtamed. In order to develop into full-scale cyclones that exhibit the usual baroclinic 
wave structure, the additional ingredient of upper-level forcing is needed. Saito (1977) 
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studied three intermediate-scale cyclones that formed during AMTEX 1975 and observed 
that these systems developed when a preexisting vortex moved off the continent into the 
East China Sea. These cyclones were strongest in the lowest layers and formed in a large- 
scale frontal zone. Development was negligible above 850 hPa. 

In this study the MCCA of the surface meridional wind also identifies an 
intermediate-scale wave disturbance over the East China Sea during northern winter. The 
characteristics of these disturbances resemble the ones studied previously, but there are 
some differences. 
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n. DATA AND METHODOLOGY 


A. DATA TYPE 

The wind data used in this study are the tropical Global Band Analyses (GBA) 
produced by the Fleet Numerical Meteorology and Oceanographic Center (FNMOC), 
Monterey, California during 1974-1988. These data are the same set used by Chang and 
Chen (1992) in their case study of the winter monsoon cold surge over the South China 
Sea. These data were produced four times daily using objective analysis on a Mercator 
grid which extends from 60°N to 40°S. Grid point spacing is 2 1/2 degrees latitude and 
longitude. The Mercator secant projection results in a change in distance between grid 
points from a minimum value of 140 km at 60°N to a maximum value of 280 km at the 
equator. The objective analysis scheme takes advantage of all the reports in the 
operational data base including surface synoptic, aircraft, pilot balloons, rawinsonde, and 
satellite data. The analysis was performed every six hours for the surface, 700, 400, 250, 
and 200 hPa levels using a successive corrections method. The six-hour persistence field 
is used as the first guess; no model forecast is used. The disadvantage in using a 
persistence field as a first guess is that the analysis suffers over data-void regions. Over 
such areas, with no new observations to update the field, persistence may last for days. 
On the other hand, the analysis is not dependent on any model dynamics. Both wind and 
temperature were analyzed and then adjusted by a set of numerical variational analysis 
equations which incorporated the dynamical constraints of the momentum equations with 
friction included in the surface layer (Lewis and Grayson 1972). 

In this study only the surface, 700 hPa, and 200 hPa 12 hourly wind data are used. 
At the surface the GBA is a particularly rich data set when compared with other 
operational surface analyses for the same time period. It includes all Navy and 
conunercial ship data as well as land-based data. The analyses are drawn closely to the 
observations as there is no model constraint. At 200 hPa the data set is enriched by 
frequent commercial aircraft reports over the tropical western Pacific. Since 1979 cloud 
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top winds derived from the Japanese Geostationary Meteorology Satellite have also been 
incorporated. The surface meridional winds are used to calculate leading MCCA modes 
that show the dominant horizontal structure of winter monsoon synoptic events. Data from 
200 hPa is used to determine associated upper-level patterns as well as calculate 
divergence which can be used as an indicator of deep tropical convection. 

The data are analyzed over a three month period that includes January, February, 
and March during the years 1975-1988. These particular three months were chosen for 
two reasons: First, we wish to minimize any effect from tropical cyclones which are 
periodically "bogussed" into the GBA data. January, February, and March typically have 
the fewest number of tropical cyclones (JTWC 1988). Also, the winter monsoon synoptic 
phenomena we wish to smdy (i.e. cold surges) occur most often during this time period. 

B. ANALYSIS METHODS 

1. Multiple-set Canonical Correlation Analysis 

The statistical method of analysis used in this study. Multiple-set Canonical 
Correlation Analysis (MCCA), is a particular type of eigentechnique designed to extract 
statistical information from datasets with large dimensions and long time periods. One of 
the earliest and most commonly used techniques of this type is the Principal Component 
Analysis (PCA) sometimes called Empirical Orthogonal Fimction Analysis (EOFA) (e.g.. 
Hotelling 1933; Lorenz 1956; Preisendorfer 1988). The purpose of PCA is to concentrate 
the internal variation (or variance) within a particular field to the leading components. The 
first principal component (PC), a function of time, is a linear combination of a given 
variable over a set of grid points that has maximum possible variance. The coefficients 
associated with each grid point used in the linear combination are known as the 
"weighting function", a function of space only. The principal component, a function of 
time only, is sometimes referred to as the "amplitude". Subsequent PC’s can be found 
under the same maximum variance requirement and under the condition that each 
weighting function is orthogonal to aU previous PC weighting functions. If a data pattern 
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is proportional to a particular PC’s weighting function, then aU other PC’s will vanish due 
to the spatial orthogonality. Therefore the weighting functions, often called EOF’s in 
meteorological studies, can be viewed as spatial patterns associated with each PC. It is 
hoped that with PCA, the leading components (those with the largest variabilities) will 
represent different physical processes. Also, component truncation of the lowest modes 
can help delete noise and errors probably associated with small variabilities. 

While PCA gives information about the variation within one data field, it offers 
no help in determining the relationship between different fields. Canonical Correlation 
Analysis (CCA) is designed to expose relationships between two fields in terms of then- 
correlation coefficients (Hotelling 1935,1936; Glahn 1968; Prohaska 1976). With CCA 
the elements of each field are linearly combined to form the canonical components. A 
pair of spatial patterns, analogous to the EOF’s in PCA, is also produced that correspond 
to each optimized correlation. Because conventional CCA only reveals maximmn 
correlations between two fields without regard to variance significance, low order mode 
truncation using principal components is often recommended before CCA is carried out 
(Barnett and Preisendorfer 1987). False relationships due to highly correlated, low 
variability noise and errors are thus avoided. 

The MCCA method is a generalization of the two field CCA to multiple fields 
(Chen et al. 1994). It seeks the optimal correlation among more than two data fields. The 
technique uses diagonalization of the product or squared product of correlation matrices 
between selected field pairs in seeking the optimal relationship between the data fields. 
The resultant canonical components approach the optimal geometrically averaged 
correlation among the data fields. Because it is unlikely for each correlation matrix to be 
diagonalized when the number of datasets is more than 2, there are non-zero cross 
component correlations called "residual correlations" for all selected pairs. If these 
residual correlations are large when compared to the canonical correlation coefficients, 
then the solutions obtained by the MCCA are doubtful. Chen and Chang (1994) 
completed a reliability case study to test the MCCA method and fotmd that it is stable 
with respect to sampling changes when the data contain significant signals of physical 
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phenomena and not stable when the data are random. 

In this study three spatial domains are used: two small core domains comprised 
of 255 grid points (15 by 17) imbedded in a large domain comprised of 1125 grid points 
(25 by 45) (Figs. 1 and 2). The first small domain covers the South China Sea and 
vicinity: 10 S-25°N, 90°E-130°E. Its choice is based on previous cold surge studies (e.g., 
Chang and Chen 1992) which show that the South China Sea is the main domain of surge 
activity. The second small domain covers the East China Sea and vicinity: 5°N-40°N, 
110E-150E. This domain is chosen to study the eastward propagating baroclinic 
disturbances that occur during the northern winter surge season. The large domain spans 
most of the western Pacific and eastern Indian oceans: 20‘'S-40°N, 70°E-180°E. The 
MCCA modes are calculated over the small domains using the north-south (v) surface 
wind component. The stmcture over the large domain is then determined by single-point 
correlation and composite analyses with respect to the basic MCC modes. Surface v was 
chosen as the variable used to compute the modes for two reasons. First, surface data is, 
in general, more reliable because of the increased availability of observations. The second 
reason is that the variations in the north-south component of the wind (v) have a stronger 
signal than the east-west component (u) over our area of study. 

Before MCCA is applied, data within the two core domains are first truncated to 
the first 18 principal components which retain 75% of the total variance for the South 
China Sea domain and 86% of the total variance for the East China Sea domain. The data 
is then processed through a 14-day high pass filter so that only phenomena with periods 
less than 14 days are retained. This filtering is necessary to ensure that the first-guess 
persistence which may mask signals of interest is removed. The MCCA is then applied 
to 12 datasets that are constructed by an overlapping sliding time window that contains 
12 consecutive v fields. The 12 datasets from 0 h (time TO) to 132 h (time Til) are 
formed as the window is sliding forward at 12 hourly intervals. Thus, the analysis looks 
for correlated variations through multiple times (Chen and Chang 1994). The squared 
product of correlation matrices between fields of adjacent times is then optimized. After 
transformation matrices from principal components to canonical components are obtained, 
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spatial patterns of the v field at each time step are computed as weighting functions that 
convert the canonical components back to the original v field. 

At any given time a mode contains both horizontal spatial and time structures, 
with the latter encompassing 12 consecutive twelve-hourly frames from day 0 to day 5.5. 
Therefore, the time variation of each mode is a space-time volume function that covers 
not only the horizontal spatial domain, but also a time domain of 5.5 days. To account 
for this time domain, a time-integrated amplitude coefficient is produced. This coefficient 
is the average amplitude of the 12 consecutive twelve-hourly structures and is plotted in 
the time domain according to the first day (day 0) time frame. 

2. Composite and Single-point Correlation Analyses 

After the MCC modes have been calculated over the small domains, the results 
can then be used to study phenomena over the large domain. One method is to use 
composites. Based on the time-integrated amphmde coefficient time series for a particular 
mode, a number of specific cases may be picked and then averaged together to give a 
composite plot. A sequence of composite plots may be produced by picking cases based 
on a given time lag or lead relative to the first plot time. Composite analysis is a zero- 
order technique; that is, we observe a mean anomaly pattern that may or may not be 
representative of any individual case. Because the cases have been selected based on a 
particular MCC mode, however, it is likely that each case is indeed similar and that the 
composite does provide a reasonable representation. 

Another method of studying the effects of the three dimensional time variability 
of the basic MCC modes on the large domain is to use single-point correlation. Here the 
amplirnde coefficient time series for a particular mode is correlated with the time series 
for each grid point in the large domain. Contours are drawn based on the grid point 
correlation values. A sequence of correlation plots based on different time lags may also 
be produced. Single-point correlation is a first-order analysis technique. The correlation 
plots produced give us some idea of how well the MCC mode is represented in the mean- 
removed data. It is expected that the correlation plot will show a similar spatial pattern 
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to the MCC mode weighting function for the same variable. 


C. DATA FILTERING 

As mentioned above, the original principal component data is processed through 
a 14-day high pass filter before MCCA is applied. Without this filtering, the leading 
MCCA modes are entirely dominated by long-period persistence that completely masks 
significant synoptic time scale events. Fig. 3 shows the weighting function of mode 1 
when MCC analysis is apphed to the East China Sea domain with no filtering. The large, 
nearly stationary feature depicted is due to a highly correlated, low-frequency background 
signal in the data field. Persistence in the original data comes from two sources. The first, 
physical persistence, is simply that component of the data that does not change from one 
time period to the next. The second source is artificial and is created when a six-hour 
persistence field is used as a first guess in the tropical GBA dataset. If physical changes 
occur in the actual field, but are not reflected in the available observation reports, 
persistence is artificially introduced. 

One method of overcoming this problem is to "time-difference" the data before 
applying MCCA. Fig. 4 shows the weighting functions of mode 1 when MCC analysis 
is applied to the East China Sea with 24-hourly differenced data; 

x’it) = xit+12 h) - x(t-12 h) 

Notice that applying time differencing acts like a high pass filter. Although this method 
is successful in eliminating the persistence problem, it is more difficult to interpret the 
weighting functions which now show first order acceleration fields rather than zero order 
velocity fields. 

The second method of overcoming this "persistence problem" is to process the 
original data through a high pass filter prior to analysis. This will eliminate the long 
period, low frequency persistent components and leave the higher frequency, synoptic- 
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scale events. Because the events examined in this study occur with a maximum period of 
about 8 to 9 days, a 14-day high pass filter was chosen. 
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Figure 1. Map showing the core South China Sea domain imbedded in the large 
western Pacific and eastern Indian ocean domain. 
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Figure 3. Weighting function for MCC mode 1, East China Sea domain, of the surface 
meridional winds (v) for 12 consecutive twelve-hourly frames from OOh (TO) to 132h 
(Til). Contour interval is 0.03 and dashed lines correspond to northerly winds when 
MCC mode 1 amplitude is positive. No prior filtering was applied to the data. 
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m. LEADING MCC MODES BASED ON SURFACE V 


A. THE SOUTH CHINA SEA DOMAIN 

The weighting functions for MCC modes 1 and 2 derived from the surface v wind 
component are shown in Figs. 5 and 6, respectively. Each figure shows the structure of 
the surface v anomaly (mean removed) for 12 consecutive twelve-hourly time frames from 
0 h (TO) to 132 h (Til) (5 1/2 days). For MCC mode 1 the average correlation between 
adjacent times is 0.82, with an average fractional variance of 21.2 % and a maximum 
residual correlation of 0.14. Recall that the smallness of the residual value compared to 
the average correlation value indicates that the solution is reliable (Chen and Chang 
1994). For MCC mode 2 the average correlation is 0.80 with an average fractional 
variance of 20.1% and a maximum residual correlation of 0.14. A lag correlation between 
modes 1 and 2 was also calculated and the results shown in Fig. 7. The lag correlation 
plot indicates that modes 1 and 2 are a pair. At a time lag of about 3.0 (1 1/2 days) the 
correlation between the two modes is approximately 0.8. We can thus conclude that 
modes 1 and 2 show basically the same phenomena and differ only by a phase shift. This 
can also be seen by noting the similarity of the weighting function patterns for modes 1 
and 2. The average fractional variance for the combined modes 1 and 2 then is 41.3%. 
The associated amplitude coefficient time series and power spectra for modes 1 and 2 are 
shown for each year in Figs. 8 through 14. 

The weighting function pattern for mode 1, time TO, is characterized by two 
large-scale areas: a positive area centered near 25'N, 103°E that covers most of the 
Indochina peninsula and a broader negative area centered north of Taiwan that covers the 
northern half of the South Qiina Sea. While the positive center descends only very 
slightly southward to near northern Thailand and weakens through time T5, the negative 
area propagates more rapidly and splits into two centers. A tight gradient over the 
Indochina peninsula is established between the positive and negative areas. This gradient 
may be due in part to the fact that the GBA siuface data used in this study are not 
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discriminant of elevation so that observations taken over mountainous regions may have 
a stronger bias compared to sea-level observations. After the broad negative area splits, 
one center moves southward into the southern South Qiina Sea while the other moves 
southeastward into the Philippine Sea. This pattern of propagation continues from times 
T6 through Til in almost the exact same fashion as just described, but with areas of 
opposite sign. This indicates that the apparent "surge" we see moves rather like a wave 
that begins in the north and dissipates in the south and not an intrusion from the north 
that advances and then retreats. Notice that the pattern at time T11 is almost exactly the 
opposite of that at time T3. This indicates a half cycle of about 4 days, so that the total 
period of the event would be about 8 days. Indeed, the average power spectrum for all 
13 years, modes 1 and 2 combined, confirms that the broad period spectral maximum is 
in a window between 6 and 12 days with a small peak at about 8 days (Fig. 15). Table 
1 is a summary of the spectral peaks for each year. While there are year to year 
variations, the peaks generally stay within the 6-12 day window and most frequently 
occur at 8 days. 

Lower order modes (modes 3, 4, etc.) were also calculated for the South Qiina 
Sea domain, but because they each contribute so little to the total variance, they are not 
included in this study. 

B. THE EAST CHINA SEA DOMAIN 

1. Modes 1 and 2 

The weighting function for MCC modes 1 and 2 derived from the surface v wind 
component over the East China Sea domain are given in Figs. 16 and 17, respectively. 
As before, the patterns show the structure of the surface v anomaly for 12 consecutive 
twelve-hourly time frames from 0 h (TO) to 132 h (Til). For mode 1 the average 
correlation between adjacent times is 0.86, with an average fractional variance of 11.3% 
and a maximum residual correlation of 0.22. For mode 2 the average correlation between 
adjacent times is 0.86, with an average fractional variance of 11.9% and a maximum 
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residual correlation of 0.20. As with the South Qiina Sea domain, MCC modes 1 and 2 
are a pair, differing only by a phase shift. This is verified upon examination of the lag 
correlation plot between the two modes (Fig. 18) which shows a correlation value of 
about 0.85 at a time lag of about 3.5 (1 3/4 days). The two modes together then accoimt 
for approximately 23.2% of the total variance of the 14-day filtered data. Again, the 
centers of the anomalies lie in a band anywhere between the latitudes 25° N and 40° N 
which indicates that MCC modes 1 and 2 here probably represent the typical observed 
wintertime barochnic waves of the western Pacific. The amplitude coefficient time series 
for modes 1 and 2 along with their power spectra for each year are shown in Figs. 19 
through 25. 

The weighting function patterns for MCC modes 1 and 2 show large anomalies 
with horizontal wavelengths on the order of 3000 to 4000 km. The anomalies have a 
northeast to southwest horizontal tilt and move from west to east with a period of about 
7-8 days. The average power spectrum for all 14 years, modes 1 and 2 combined, is 
strongest between 6 and 12 days with a significant peak at about 8 days (Fig. 26). Table 
2 is a summary of the spectral peaks for each year. Again there are year to year 
variations, but the peaks generally fall within the 6-12 day window and most frequently 
at about 8 days. 

The horizontal spatial and time scales represented in modes 1 and 2 for the East 
China Sea domain look very similar to those in the South China Sea domain. Therefore, 
although calculated over different spatial domains, both pairs of modes probably are 
manifestations of the same event. To confirm this, coherence and phase calculations were 
carried out between the two pairs’ amplitude coefficient time series (Fig. 27). Significant 
coherence values that range between about 0.65 and 0.80 are found for the 6 to 12 day 
period. The transfer function phase over this same period is nearly out of phase 
throughout the frequency domain. One may conclude then that the two mode pairs are 
essentially the same, except that the patterns isolated are out of phase. Therefore, only 
MCC mode 1 for the South China Sea domain will be studied. 
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Table 2: Significant power spectrum peaks for the East China Sea modes 1 and 2 



































2. Modes 3 and 4 

Unlike the South China Sea domain, MCC modes 3 and 4 for the East China Sea 
domain contain a significant percentage of the total variance. The weighting functions are 
given in Figs. 28 and 29, respectively. The associated amplitude coefficient time series 
and power spectrum plots for each year are given in Figs. 30 through 36. Because the 
phenomenon isolated in modes 3 and 4 is well correlated over a shorter time period, a 
6-ffame (2.5 days) overlapping sliding window was used. For mode 3 the average 
correlation between adjacent times is 0.79, with an average fractional variance of 9.8% 
and a maximum residual correlation of 0.14. For mode 4 the average correlation between 
adjacent times is 0.73, with an average fractional variance of 6.2% and a maximum 
residual correlation of 0.27. Examination of the lag correlation plots given in Fig. 37 
reveals that modes 3 and 4, like modes 1 and 2, are a pair with a correlation value of 
about 0.8 at about 2.0 time intervals (1 day). As a pair, modes 3 and 4 accoimt for 16% 
of the total variance of the 14-day filtered data. The plots also show that the two mode 
pairs isolate separate events; correlation values between the two pairs are low (between 
modes 1 and 3 for example). 

The weighting function patterns for MCC modes 3 and 4 are characterized by 
elongated anomalies that tilt from northeast to southwest and have horizontal wavelengths 
on the order of 1500 to 2(XX) km. The centers of the anomalies lie in a band of latitudes 
between 20 N and 40 N and move to the east-northeast with a period of about 3 days. In 
general, the weighting functions for this pair of modes looks somewhat similar to those 
of the modes 1 and 2 in that they all show eastward propagating wave patterns. The 
primary difference appears to be in the spatial and temporal scales both of which are 
smaller for the mode 3 and 4 pair. Table 3 provides a summary of the significant spectral 
peaks for each year. Every year has a peak within a 2-4 day window, with most occurring 
right at 3 days. A few years also have peaks at longer periods. Comparison of the average 
power spectrum plot for modes 3 and 4 (Fig. 38) with that of modes 1 and 2 (Fig. 26) 
reveals that the longer time scale events of modes 1 and 2 and the shorter time scale 
events of modes 3 and 4 do not overlap. There is a spectral gap between 4 and 6 days. 
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L= peak is large relative to other years 
M= peak is moderate relative to other years 
S= peak is small relative to other years 

































































Figure 7. Lag-correlation plots for MCC modes 1 and 2, South Qiina Sea domain. 
The two modes correlated in each plot are labeled in the upper left corners. Each 
interval on the x-axis represents 2.0 days (4*12h). 
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Figure 8. Time series of the time-integrated amplitude coefficients of MCC modes 1 
and 2, South China Sea domain, for the winters 1975 and 1976. The power spectrum 
for each time series is shown at the right with the 3, 6, 9, and 12 day periods labeled. 
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Figure 10. Same as Fig. 8 except for winters 1980 and 1981. 
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Figure 15. Average power spectrum for 13 winters, MCC modes 1 and 2 combined. 
South China Sea domain, with the 6, 9, and 12 day periods labeled. 
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Figure 16. Weighting function tor MCC mode 1, East China Sea domain, of *e 
surface meridional wind (v) for 12 consecutive twelve-hourly frames from OOh (TO) to 
132h (Til) Contour interval is 0.03 and dashed lines correspond to norther y w 
when MCC mode 1 amplitude is positive. Data was passed through a 14-day high pass 

filter prior to analysis. 
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Figure 18. Lag correlation plots for MCC modes 1 and 2, East China Sea domain. 
The two modes correlated in each plot are labeled in the upper left corners. Each 
interval on the x-axis represents 2.0 days (4*12h). 
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Figure 20. 
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Figure 25, Same as Fig. 19 except for winters 1987 and 1988. 
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Figure 26. Average power spectrum for 14 winters, MCC modes 1 and 2 combined. 
East China Sea domain, with the 6, 9, and 12 day periods labeled. 
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Figure 27. Coherence and transfer function phase plots for the amplitude coefficient 
time series of MCC mode 1, South China Sea domain, and MCC mode 1, East China 
Sea domain. 3, 6, and 12 day periods are labeled. 






Figure 28. Weighting function for MCC mode 3, East China Sea domain, with the 
surface meridional wind (v) for 6 consecutive twelve-homly frames from OOh (TO) to 
60h (T5). Contour interval is 0.03 and dashed lines correspond to northerly winds 
when MCC mode 3 is positive. Data was passed through a 14-day high pass filter 
prior to analysis. 
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Figure 32. Same as Fig. 30 except for winters 1979 and 1980. 
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Figure 33. Same as Fig. 30 except for winters 1981 and 1982. 
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Figure 34. Same as Fig. 30 except for winters 1983 and 1984. 
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Figure 37. Lag correlation plots for MCC modes 1,2,3, and 4, East China Sea domain. 
The two modes correlated in each plot are labeled in the upper left comers. Each 
interval on the x-axis represents 2.0 days (4*12h). 
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Fiiure 38, Average power spectrum for 14 winters. MCC modes 3 and 4 combined, 
East Oiina Sea domain, with the 2 and 4 day periods labeled. 
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rV. STRUCTURE OF MOTIONS BASED ON CORRELATION AND 

COMPOSITE ANALYSES 

A. SOUTH CHINA SEA MCC MODE 1 

1. Single-point Correlation Analysis 

The single-point correlation analysis plots of the surface rtieridional wind (v) with 
MCC mode 1, South China Sea domain, is given in Fig. 39. As expected, the surface v 
correlation patterns over the core domain look very similar to the MCC mode 1 weighting 
function. Recall that the MCC analysis was performed using surface v data. Over the 5 
1/2 day period we again see a "surge cell" move southward into the South China Sea and 
dissipate by T7 followed by a similar cell of opposite sign. The tight gradient over the 
Indochina peninsula found in the weighting function patterns is also present here. The 
correlation pattern over the large domain contains large-scale anomalies with wavelengths 
on the order of 3000-4000 km that have northeast to southwest horizontal tilts and move 
from east to west with a period of about 8 days. As mentioned previously, these 
anomalies probably represent typical wintertime baroclinic waves. As further confirmation 
that MCC mode 1 for the South China Sea domain is essentially the same as MCC mode 
1 for the East China Sea domain, surface v data is correlated with the latter and plotted 
in Fig. 40. Both plots contain the same basic features. 

The single-point correlation analysis plots of the 700 hPa meridional wind with 
MCC mode 1, South China Sea domain, is given in Fig. 41. Here the correlation values 
are slightly smaller ranging from -0.3 to +0.3 compared to -0.4 to +0.5 for the surface 
case. At 700 hPa the core domain also shows a negative phase anomaly propagate 
southward into the South China Sea through T5 followed by a positive one through Til. 
Unlike the surface, however, the anomalies do not penetrate as far south to the equator. 
The primary feature in the large domain at 700 hPa is the 8-day, 3000-4000 km wave that 
propagates eastward from the Asian continent. The 200 hPa meridional wind correlation 
with MCC mode 1 also shows a longwave pattern with similar characteristics indicating 
that the phenomenon has large vertical extent (Fig. 42). In this case correlation values 
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range from -0.3 to +0.3 and there is no significant surge signal over the South China Sea. 


2. Composite Analysis 

a. Strong Cases 

The 13 strongest cases of cold surge over the 14 year period from 1975- 
1988 were selected based on the highest peaks of the time-integrated amplimde coefficient 
time series for MCC mode 1 shown in the previous chapter (Figs. 8-14). Surface vector 
plots over the large domain based on composites of these 13 strong cases are given in 
Fig. 43. Each of the 6 frames is separated by 24 hours (vice 12 hours as in previous 
charts). The sequence shows a strengthening of northeasterly winds over the South China 
Sea from TO to T8 and then a slight weakening at TIO. The subtropical ridge pushes 
eastward from mainland China across the East China Sea at TO, breaks off as a closed 
anticyclone or "bubble high" by T2, and then continues its eastward propagation through 
TIO. Notice that the strengthening of the northeasterly winds is mainly confined to the 
South China Sea region. 

In order to see more of the detail hidden in this sequence of charts, the 
composites are calculated for 12 twelve-hourly frames with the background mean 
removed so that just cold surge anomalies are depicted (Fig. 44). The most prominent 
feature at TO is the large cyclonic circulation anomaly centered over China producing 
strong southwesterly flow over the northern South China Sea and northwesterly flow over 
Indochina. A weak ridge sits just north of the cyclone near Mongolia. Weak anticyclonic 
anomalies are present over the Tibetan plateau and just east of Japan. Through T5 the 
large cyclonic anomaly over China becomes elongated into a trough and moves eastward 
over the western Pacific. Behind it, the weak Mongolian ridge anomaly strengthens and 
becomes a closed anticyclone. The weak Tibetan plateau cyclone loses its identity while 
the anticyclone east of Japan moves further eastward. The strong southwesterly flow 
anomaly over the South China Sea weakens considerably and is beginning to be replaced 
by northeasterlies descending from the north. During times T6 through T11 we observe 
almost exactly the opposite sequence of events. The anticyclone centered over central 
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China at T6 becomes elongated and moves eastward to just south of Japan by Til. It is 
replaced by a broad-based trough that moves in over China from the northwest. 
Northeasterly flow strengthens over the South China Sea through T8 and then weakens 
through Til. Although the magnitude is small, there may also be some evidence of cross- 
equatorial influence of the surge at times TIO and Til. As the surge weakens over the 
South China Sea, there is an increase in northeasterlies off the northwest coast of 
Australia. 

The 700 hPa vector composite plots for the same 13 strong cases are shown 
in Fig. 45. Again the background mean has been removed to reveal detail in the 
anomalies. Over the South China Sea core region we again see southwesterly flow 
weaken and shift to northeasterly flow, however the transition occurs about 24 hours 
earlier and the overall magnitude is smaller than at the surface. In the 700 hPa large 
domain a large cyclonic circulation over China at TO dissipates by T3 as a ridge moves 
in from the west and takes its place. The ridge becomes a weak semi-closed circulation 
at T5. Meanwhile, a large anticyclone centered over Japan moves eastward into the 
western Pacific and is replaced by a cyclonic circulation at T7. The cyclone follows the 
anticyclone out to sea through Til ahead of another ridge that develops over the Yellow 
Sea. As seen at the surface, there again appears to be possible evidence of cross- 
equatorial surge influence during time TIO and Til when northwesterly flow across the 
northwest coast of Austraha increases in response to a developing cyclone over the Timor 
Sea. 

The 200 hPa anomaly vector composite plots for the strong cases depict 
the expected longwave pattern propagation across the subtropical latitudes (Fig. 46). Over 
the core domain weak anomalous flow from the northeast begms about 48 hours earlier 
than at the surface and shifts to southwesterly flow at T7. One particularly interesting 
feature of this 200 hPa composite sequence are the circulations north of the Bay of 
Bengal. A weak anticyclone centered near Bangladesh at T1 strengthens through T4 and 
then becomes a ridge that drifts eastward across southern mainland <2hina and then 
northeastward along the China coast to the East China Sea by Til. Meanwhile, over 
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Bangladesh, a cyclonic circulation is firmly established by T8 and remains relatively 
stationary while strengthening through Til. 

Whether these circulations are a subtropical response to the cold surge or 
simply a reflection of the midlatitude triggering mechanisms is unknown and warrants 
further study. Chang and Lau (1982) studied 200 hPa winds of foxu winters and found a 
relationship between the acceleration of midlatitude jet streaks and the cold surge, but did 
not observe these circulations. The reason for this is that the circulations are not 
detectable when masked by the mean flow. Fig. 47 shows the same 200 hPa composites 
with the mean flow included. Here the time difference between each frame is 24 hours. 
Propagation of jet streaks across the midlatitudes is clearly shown, but there is no 
evidence of the circulations. There is also no cold surge signal detectable over the South 
China Sea as was the case with the anomaly plots. 

b. Medium Cases 

In order to determine whether the features of the cold singe described in 
the previous section are common not only to a few strong events over a 14 year period, 
but also to cold surges in general, the time-integrated amplitude coefficient time series 
for MCC mode 1 was again used as a basis to select 54 "medium-strength" cases over the 
same period. Vector anomaly composite plots were again produced for the surface, 700 
hPa, and 200 hPa levels (Figs. 48, 49, 50, respectively). Most of the characteristics of the 
surge found in the strong cases are also present in the medium cases with a few subtle 
differences. At the surface the magnitude of the northeast and southwest wind anomalies 
over the South China Sea are slightly smaller in the medium case composites. Also, 
possible cross-equatorial surge influence is manifested differently with the medium cases. 
If we focus on an area just north of 20°S between 100°E and 110”E, we notice slight 
southerly flow from TO through T5 transition to northerly flow through Til. Although 
the magnitude of these anomalies is small, the manner and period of change exactly 
match those of the cold surge. 

At 7(X} hPa the patterns for the medium cases are essentially the same as 
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those for the strong cases. Again, the difference between the two is magnitude, with the 
strong case anomalies being slightly larger. Also, the northeast and southwest anomalies 
over the South China Sea do not reach quite as far south in the medium cases, so cross- 
equatorial influence does not appear to occur. 

At 2(X) hPa the differences between the strong and medirmi cases are again 
subtle. Although the anomalies in the medium cases are smaller in magnitude, they still 
show the same midlatitude longwave propagation pattern. The anomalous circulations 
centered over Bangladesh found in the strong cases are also present here. 

c. Divergence 

As mentioned in the introduction, cold surge events may have a significant 
impact on tropical convection. To investigate this, 200 hPa divergence fields were 
cailculated over the large domain using the u and v anomaly fields. The following 
calculation was performed at each grid point: 




u(,i+l,j) - uii-l,j) v(i, 7 + 1 ) - v(i, 7 -l) 


Ax 


Ay 


The resulting divergence fields were then spatially averaged and composited based on 
the time-integrated amplitude coefficient time series for MCC mode 1 exactly as before. 

The composite sequence based on the 13 strongest cases is given in Fig. 
51. The sequence shows three noticeable features. The first is in the area defined by 0' 
to 20° S and 70°E to 90°E. This region contains a convergence anomaly from TO to T5 
that is strongest at T2. At T6 a positive divergence anomaly becomes established and 
persists through TIO. The second notable feature is found in the vicinity of Borneo. Weak 
200 hPa convergence over Borneo between TO and T7 gives way to strong divergence 
through T9. The divergence center is particularly strong at T9. In both areas the increases 
in 200 hPa divergence is evident in the second half of the sequence. This coincides with 
the time of maximum northeast surge anomaly over the South China Sea described in the 
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previous section and is consistent with the findings of Lau and Chang (1987). The 
difference in Lau and Chang’s study is that they used velocity potential fields which 
smooth out any synoptic scale divergence. Our results show anomalies with a much 
smaller zonal scale. The third interesting feature is a positive divergence anomaly 
positioned along the equator near ISO'E at time TO. This positive anomaly persists 
through the entire 12 frame sequence and slowly moves eastward to near 175’E by Til. 
This feature might be associated with an equatorial Kelvin wave. 

A composite sequence of 200 hPa divergence fields based on the 54 
medium-strength cold surge cases was also produced (Fig. 52). The most striking 
difference between this and the strong cases sequence is the lack of significant magnitude 
anomalies. In fact, the maximum values shown are probably well within the range of error 
inherent in kinematic divergence calculations. We might conclude then that the medium- 
strength cold surge does not produce enough of an impact on 200 hPa divergence to be 
detectable in the data. 

B. EAST CHINA SEA MCC MODE 3 

1. Single-point Correlation Analysis 

The single-point correlation analysis plots of the surface meridional wind (v) with 
MCC mode 3, East China Sea domain, is given in Fig. 53. Again, as expected, the surface 
V correlation patterns look very similar to the MCC mode 3 weighting function. 
Correlation values range from -0.4 to +0.4. Over the 2 1/2 day sequence, elongated 
anomalies that tilt from northeast to southwest with wavelengths on the order of 1500 to 
2000 km move generally eastward with periods of about 3 days. The anomalies appear 
to develop over the Asian continent and then move out into the East China Sea before 
tracking more northeasterly to the north-central western Pacific. 

The single-point correlation analysis plots of the 700 hPa meridional wind with 
MCC mode 3 is given in Fig. 54. Again correlation values range from -0.4 to +0.4. The 
pattern here looks exactly as the surface case except the anomalies are more round and 
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do not extend as far southward. The 200 hPa meridional wind correlation with MCC 
mode 3 also shows the same characteristic anomalies, but here the centers are displaced 
to the north along 40*N (Fig. 55). Correlation values at 200 hPa also range from -0.4 to 
+0.4. 


2. Composite Analysis 

The time-integrated amplitude coefficient time series for MCC mode 3 over 14 
winters from 1975 to 1988 (Figs. 30-36) was used as a basis to select cases for composite 
analysis. 27 cases were selected and composited to yield 6 twelve-hourly surface anomaly 
vector plots shown in Fig. 56. The sequence is characterized by a trough at TO centered 
over the Sea of Japan that propagates eastward over the central western Pacific by T5. 
Behind it, a ridge anomaly located along the China coast at TO becomes a closed 
anticyclone near Korea at T2 and moves almost due east to near 150°E by T5. One 
particularly interesting feature in this sequence is the reversal of wind anomalies over the 
South China Sea region from northerlies between TO and T2 to southerlies from T3 
through T5. This indicates that there is also some effect, albeit small, over the South 
China Sea area associated with smaller, shorter time scale events in the midlatitudes. This 
implies that the strong cold surge over the South China Sea may be expected when MCC 
modes 1 and 3 are in phase, and weak surges may be expected if they are out of phase. 

The 700 hPa anomaly vector composite plots for the same 27 cases are shown in 
Fig. 57. The same features seen in the surface plots are also present here except they are 
displaced slightly to the west indicating they have a backward tilt with height, which is 
characteristic of barochnic systems. There is no apparent signal over the South China Sea 
region at 700 hPa. 

At 200 hPa, the anomalies again have the same spatial and temporal scale as the 
lower levels with a small displacement to the west (Fig. 58). The reversal in wind 
anomalies over the South China Sea region seen in the surface plots is also present here. 
Note that southwest anomalies at TO appear to be a weak southwestern extension of the 
midlatitude jet streak centered over Japan. As the jet streak moves east of Japan through 
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T2, the magnitude of the anomalies decreases to near zero. This phenomenon is also 
evident in the 200 hPa MCC mode 1 composite sequences for both the strong and 
medium-strength cases although it is not as obvious. 

The fact that these shorter wavelength, shorter time scale anomalies are present 
at 700 hPa and 200 hPa and have westward vertical tilt indicates they are baroclinic with 
large vertical extent. Recall that the intermediate-scale cyclones observed by Nitta and 
Yamamoto (1973) discussed in Chapter I, were of the same horizontal spatial and 
temporal scales as these here, but are described as being shallow with no sigmficant 
development above about 850 hPa. It is possible that the shorter scale disturbances 
represented here begin as weak shallow depressions over the east Asian continent and 
then gain upper level support as they move out over the ocean. This scenario is consistent 
with the findings of Saito (1977) mentioned earlier. 
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Figure 39. Single-point correlation patterns of surface v with MCC mode 1, South 
China Sea domain, at 12 consecutive 12-hourly time lags. Contour interval is 0.1 and 
negative values are dashed. 
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Figure 48. Surface wind vector plots based on composites of 54 medium-strength cold 
surge cases at 12 consecutive 12-hourly time lags with the background mean removed. 
Vector scale is shown at the top of each frame. 
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Figure 56. Surface wind vector plots based on composites of 27 cases at 6 consecutive 
12-hourly time lags with the background mean removed. Vector scale is given at the 
top of each frame. 
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V. SUMMARY AND CONCLUSIONS 


In this study we examine the structures of motions over a large domain that 
includes the western Pacific and eastern Indian oceans dtuing 14 winters from 1975-1988 
using Navy tropical Global Band Analysis (GBA) data. Multiple-set Canonical Correlation 
Analysis (MCCA) is applied to 12 datasets in two small domains that are imbedded in 
a large domain that includes the western Pacific and eastern Indian oceans. The structure 
of motions in the large domain is then determined through the use of single-point 
correlation analysis and composite analysis based on the basic MCC modes. The 12 
datasets are constructed by a forward sliding time window which consists of 12 twice- 
daily time frames of the surface v field, through 3 months of the year (January, February, 
March). The resulting MCC modes are ranked by the geometrically-averaged correlation 
between surface v fields of adjacent times within the 12 time frames, or 5.5 days. 

MCC modes 1 and 2 calculated over the small South China Sea and East China 
Sea domains reveal the basic stmcture of the Asian winter monsoon cold surge. The surge 
moves across the South China Sea with a positive anomaly moving from north to south 
followed by a negative anomaly with an average period of 8 days. Previous studies of the 
cold surge (e.g., Chang and Chen 1992) did not address the character of the "break" in 
the cold surge. Our result indicates that the break does not occur as a withdrawal of the 
northeasterly wind northward, but rather as an equatorward propagation of the negative 
phase. Meanwhile, large-scale disturbances with wavelengths on the order of 3000 to 4000 
km propagate eastward from Asia across the East China Sea also with an average period 
of 8 days. Coherence and phase calculations between the amplitude coefficient time series’ 
of MCC mode 1 for the South China Sea domain and MCC mode 1 for the East China 
Sea domain confirm that they are both manifestations of the same event (i.e. the cold 
surge). In addition, the time series’ are also used to construct single-point correlation 
analyses over the large domain. Both MCC mode 1 time series’ yield similar correlation 
patterns. In summary, the surge can be viewed as a single mode over a rather broad 
region with an average period of 8 days that includes eastward propagating longwaves 
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over the East China Sea and equatorward propagating surge waves across the South China 
Sea. 

Other features evident during cold surge events are depicted in composite plots 
over the large domain. Cases are selected based on the magnitude of the peaks of the 
time-integrated amplitude coefficient time series for MCC mode 1. Over the 14-year 
period 13 cases are categorized as "strong" and 54 cases are categorized as "medium- 
strength". Goss-equatorial surge influence is seen in both the strong and medium-strength 
cold surge surface composites. For the strong cases, an increase in northeasterly winds 
off the northwest coast of Australia is exhibited during the maximum northeasterly surge. 
For the medium-strength cases, a phase shift in anomalous winds with a period that 
exactly matches that of the cold surge is observed in the area between 100°E and llO'E 
along 20'S. The cross-equatorial influence of cold surges has been recognized by other 
investigators (e.g., Lau 1982; Lau and Chang 1987). At 200 hPa, composite plots reveal 
a large-scale anomalous anticyclone positioned north of the Bay of Bengal that is 
gradually replaced by an anomalous cyclone throughout the surge sequence. This feature, 
which was not reported in previous cold surge studies, exists for both the strong and 
medium-strength composites. As mentioned previously, whether the circulations are a 
cause or effect of the cold surge remains unknown. Perhaps reanalysis over a larger 
domain that includes higher latitudes would provide some insight. 

Several investigators have found that there is an increase in the level of convective 
activity over the general area of the southern South China Sea during monsoon surges 
(e.g., Houze et al. 1981; Chang et al. 1979). In order to determine whether or not 
evidence of this exists in the GBA data, divergence fields are calculated over a tropical 
band between 20'S and lO'N using 200 hPa u and v components and then composited 
based on strong and medium-strength amplitude coefficients. Results show a noticeable 
increase in divergence over Borneo for the strong cases, but not for the medium-strength 
cases. This indicates that only during strong cold surges is the 200 hPa divergence strong 
enough to be detectable by our kinematic calculation method. The divergence field 
composites also show a positive anomaly located along the equator near New Guinea at 
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the time of the surge onset that slowly moves eastward toward 180° over the 5 1/2 day 
sequence. This anomaly might be associated with an equatorial Kelvin wave. 

Over the East China Sea MCC modes 3 and 4 reveal eastward propagating 
anomalies with wavelengths on the order of 1500 to 2000 km with an average period of 
about 3 days. Composite plots of 27 cases selected from the 14 winter seasons show a 
sequence of shortwaves propagating across the midlatitudes that look similar to the MCC 
mode 3 weighting function patterns. A reversal in wind anomalies from northerly to 
southerly with a period of 3 days is present across the South China Sea. This "mim surge 
phenomenon present in MCC modes 3 and 4 will provide a contribution to the primary 
surge seen in MCC modes 1 and 2 when the two pair are in phase. Nitta and Yamamoto 
(1973) also found cyclones that propagate across the East China Sea with similar 
wavelengths and periods, however they were shallow with little development above 850 
hPa. Our results show that these disturbances are in fact deep. The 200 hPa composite 
sequence for the same 27 cases show 3-day shortwaves moving off the Asian contment 
that are about a quarter cycle out of phase with the surface disturbances. 

Future investigations of the East Asian cold surge may be conducted with similar 
statistical techniques used in this study, but with new, improved datasets. The GBA data 
used here, although it covers a relatively long period, suffers from an analysis scheme that 
uses persistence as a first guess. Over the data-sparse regions typical in the tropics, this 
artificial persistence may last for days. Major numerical weather prediction centers in the 
United States and Europe are now producing reanalyzed data that goes back to the late 
1970’s. These reanalyses are four-dimensional assimilated data that use model forecasts 
as a first guess rather than persistence fields. In these reanalyzed data all variables at all 
standard levels are available for a detaUed analysis. This may offer more insight into 
some of the features associated with cold surges that were examined in this study. 
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